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Abstract 

We prove exponential localization for the Schrodinger operator with a Poisson random potential at the bottom 
of the spectrum in any dimension. We also prove exponential localization in a prescribed interval for all large 
Poisson densities. In addition, we obtain dynamical localization and finite multiplicity of the eigenvalues. 

Resume 

On demontre localization exponentielle pour l'operateur de Schrodinger avec un potentiel aleatoire de Poisson, 
pour les basses energies et en toute dimension. On demontre aussi localization exponentielle dans un intervalle 
d'energies donne et a grande densite. On obtient de plus localisation dynamique et le fait que la multiplicite des 
valeurs propres est finie. 



1. Results 

The Poisson Hamiltonian is the random Schrodinger operator on L 2 (R d ) given by 

H x = -A + V x , with Vx(x) = ^2u(x-C), (1) 

where the single-site potential u is a nonnegative C 1 function on R d with compact support-without loss 
of generality we take u(x) < u(0) = 1-and Vx is a Poisson random potential, that is, X is a Poisson 
process on M. d with density g > 0. Thus the configuration X is a random countable subset of M. d , and, 
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letting Nx(A) denote the number of points of X in the Borel set A C I 11 , each Nx(A) is a Poisson 
random variable with mean g\A\ (i.e., F g {N x (A) = k} = (g\A\) k (k\)- 1 e'^ A ^ for k = 0,1,2,...), and 
the random variables {Nx{Aj)}™ =l are independent for disjoint Borel sets {A,}" =1 . We will denote by 
(X,F g ) the underlying probability space for the Poisson process with density g. 

Note that H x is an ergodic (with respect to translations in E d ) random self-adjoint operator. It follows 
that the spectrum of H x is the same for P e -a.e. X , as well as the decomposition of the spectrum into pure 
point, absolutely continuous, and singular continuous spectra. For u as above we actually get <j(H x ) = 
[0,+oo[ for P e -a.e. X [9]. 

We prove exponential localization for Poisson Hamiltonians at the bottom of the spectrum. By xb we 
denote the characteristic function of the set Bcl d , with \x denoting the characteristic function of the 
cube of side 1 centered at x G M d . We write (x) = yl + \x\ 2 , T(x) = {x) v for some fixed v > | . 

Theorem 1 Given g > 0, there exists Eq = Eo(g) > and m — m(p) > 7 such that for ¥ g -a.e. 
X the following holds : the operator Hx has pure point spectrum in [0, Eq] with exponentially localized 
eigenf unctions with rate of decay m, i.e., if '(f) is an eigenfunction of H x with eigenvalue E G [0, Eq], there 
is a constant < oo such that 

WxAW < C e- m|a:| for all x G R d . (2) 

Moreover there exist constants t > 1, s <=]0, 1[, and C < oo, such that for eigenfunctions ip,(f> (possibly 
equal) with eigenvalue E G [0, Eq] we have 

IIXsV'll ||X»0II < CWT- 1 ^ HT-Vll e^ T e-\ x ~y\ 3 for all x, y E Z d . (3) 

In particular, the eigenvalues of H x in [0, Eq] have finite multiplicity, and Hx exhibits dynamical local- 
ization in [0, Eq], that is, for any p > we have 

sup\\(xye-* tH x X[0iEa] (Hx)xo\\i < ™- (4) 
t 

For Poisson random potentials the density g is a measure of the amount of disorder in the medium. 
The next theorem gives localization at high disorder. 

Theorem 2 Given E n > 0, there exists g n > such that for g > go the conclusions of Theorem 1 hold 
in the interval [0, Eq]. 

While Poisson Hamiltonians are the most natural random Schrodinger operators in the continuum (the 
distribution of impurities in a material being naturally modeled by a Poisson process), a mathematical 
proof of the existence of localization has been a long-standing open problem. Localization has been known 
only in one dimension [11]. A Poissonian model, which incorporates random intensities with bounded 
densities and requires single-site potentials that do not decay too slowly at infinity, was considered in [3] . 

In the multi-dimensional case, localization in the continuum had been proved for Anderson-type Hamil- 
tonians with random intensities with bounded densities, e.g. [3], and for an R d -crgodic Schrodinger oper- 
ator with a Gaussian random potential [5] ; in both cases there is an "a priori" Wegner estimate obtained 
by averaging with bounded densities. But recently Bourgain and Kenig [2] proved localization for the 
Bernoulli- Anderson Hamiltonian, with the Wegner estimate being proven in a multiscale analysis. 

To prove Theorems 1 and 2 we exploit the new ideas introduced by Bourgain and Kenig [1,2]. In 
particular, the control of the resonances (the Wegner estimate) is achieved by a multiscale analysis 
using "free sites" and a new quantitative version of unique continuation which gives a lower bound on 
eigenfunctions. 

The control on the eigenfunction correlations given in (3) was introduced in [8]. That (3) implies 
dynamical localization is rather immediate. As for the finite multiplicity property, it follows by estimating 
\\XxX{E}{Hx)\\l\\XvX{E}(Hx)\\l from (3) and summing over x G Z d . 
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In the next section we outline the main ideas in the proof of Theorem 1 ; detailed arguments will be 
given in [6]. Theorem 2 is proved in a similar way, although the proof requires some modifications. 



2. The main ideas 

Given a cube A in R d , we let Hx.a = — Aa + Vx,a be the restriction of H\ to A with Dirichlet 
boundary condition. We consider the finite volume resolvent R\(E) = (Ha — E)^ 1 (we will omit X from 
the notation). By c, c', . • • we denote positive constants (not necessarly the same) independent of g, L, . . . . 

The multiscale analysis requires an a priori probabilistic estimate on the resolvent Ra (E) for all E G 
[0, E ] at a fixed, but sufficiently large, initial scale L , where A is a cube of side L centered at, say, 
xq. To obtain this initial estimate for Theorem 1, we divide the cube A into non-overlapping cubes 
A(j) of side £$ ~ (g~ 1 logL ) 3 centered at points j e x a + £ Z d . We consider configurations such 
that 1 < N(A(j)) < g£ d for all A(j), an event with high probability, more precisely, with probability 

— 1 — (^) ^ P ~ ^ ~ Lq P+(1 , where we can arrange for p large as desired. 

For such configurations, we pick one Q e A(j), and split the potential as Va = + V A \ with 
V^ 1 ' = J2'j u ( x ~ 0)) where J^'j denotes the sum over sites j € xq + 2£oZ d only; as a consequence the 
u(x — Q) in the sum arc non-overlapping. We have < < 1 and V^ 2 ' > 0. In order to estimate 
||-Ra(-E)||, it is convenient to use the operator T\(E), defined by 

Ta(E) = (H {2) + 1)^(1 + E-V [ A 1) )(H {2) + with H {2) = -A A + V (2) > 0. (5) 

Proceeding as in [2, Section 4], suppose ||Pa(-E)| > 1 — E a with E small. Then there is g € L 2 (M d ), with 
1 - V^b < \\g\\ < 1 and, using V^ 2) > , ||V.g|| < 2£ 3 , such that for each a € A we have 

0< (raV^g^) <cE$(\a\ + l), (6) 

where r a denotes translation by a and and the estimate is uniform in L (cf. [2, Eqs. (4.7), (4.8), and 
(4.10)]). On the other hand, taking K = 104), and recalling the definition of V A \ we get (cf. [2, Eqs. 
(4.12) and (4.15)]) 

/ r a (vi 1) )da>c X A L withoO. (7) 

J[-K,K] d 

Combining (6), (7), and the lower bound on \\g\\, we get c(l — ^/Eq) 2 < d ' E^ K d+1 , which leads to a 
contradiction for E n w £ ~( 4 ( d+1 )+) and L large. 

We may thus conclude that if g is fixed, p > is given, E G [0, E ] with E ~ (g~ x log L a )~( a ^ , and 
Lo is sufficiently large, then, with probability > 1 — Lq P , we have ||i?A(-E)|| ^ Eq 1 and \\xxRA(E)xy\\ ^ 
e -cL f or x ^ y g yY with | a; — y\ > jft. Moreover, it is clear that if = X^ey u ( x ~~ C)i the results are 
still valid if we replace by J2<;eY t( u ( x ~ with arbitrary € [0, 1]. We now declare all boxes A(j) 

with j £ x + 2£ Z l d (and hence do not contribute to V^) to be free boxes. Moreover, inside the free 
boxes we use the representation of the Poisson process X by a thinned Poisson process (e.g., [10]), that 
is, by a Poisson process Y with density 2g in such a way that to each Poisson point £ G Y is attached a 
Bernoulli random variable = or 1 with equal probability, and the single-site contribution to the 

potential is given by e^u(x — £). Note that any site £ <G Y" in a free box is a free site in the sense of [2]. 

The multiscale analysis now proceeds by induction. If A is a box of size L, we divide it into non- 
overlapping cubes A(w) of side ss e~ L centered at points w € e~ L Z d . With probability > 1 — L~ p , p 
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large, we require N(A) < gL d and all N(A(w)) < 1. We introduce an equivalence relation on Poisson 
configurations X\ in A ; X\ is the collection of Poisson configurations in A that cannot be distinguished 
from X\ by the counting functions N(A(w)). The crucial observation is that if we change a Poisson 
configuration to another one in the same equivalence class, then the eigenvalues of Ha in a fixed interval 
do not move by more than < e~ L . We may thus consider only the case when the Poisson points in A 
are in the lattice e n Z , since the desired results will then hold for the whole equivalence class. This 
reduction allows the use of the results in [2] , using equivalence classes of Poisson configurations instead 
of fixed Bernoulli configurations. Inside the free boxes equivalence classes are defined as above but for 
the Poisson points in Y. Since we have a finite number of equivalence classes of configurations inside a 
free box, we fix the points of the Poisson process Y in the free boxes, and conduct the analysis of [2, 
Lemmas 5.1], "tuning" the free parameters t% to = or 1 to obtain "good" configurations, with a 
probability estimated by Sperner's Lemma using [2, Lemma 3.1]. As in [2], we get the following result 
(cf. [2, Proposition A]), where A^ denotes a cube of side L. 

Proposition 1 Given g > 0, there exists E n = E (g) > and L n — L (g) < oo, such that if X\ L (E) 
denotes the Poisson configurations for which 

\\Ra l (E)\\ < e Ll ~ and \\ Xx Ra l (E) Xy || < e~ cL for \x - y\ > A ( 8 ) 

then for all L > L and all E e [0, E n ] we have 

F e {X AL (E)}>l--^-. (9) 

Proposition 1 provides a single-energy multiscale analysis. The weak probability estimate in (9) does 
not allow for an energy- interval multiscale analysis as in [4,7]. The first part of Theorem 1, namely 
exponential localization, requires the energy elimination scheme given in [2, Section 7]. To obtain the 
decay of the eigenfunction correlations given in (3) we add ideas from [8] . 
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